Quantum fragile matter: mechanical excitations of a Reggeon ion chain 
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This paper proposes to study quantum fragile materials with small linear elasticity and a strong response 
to zero-point fluctuations. As a first model, we consider a non-unitary (but !P7~- symmetric) massive quantum 
chain with a Reggeon-type cubic nonlinearity. At the critical point, the chain supports neither the ordinary 
quantum phonons of a Luttinger liquid, nor the supersonic solitons that arise in classical fragile critical points 
in the absence of fluctuations. Quantum fluctuations, approximately captured within a one-loop renormalization 
group, give rise to mechanical excitations with a nonlinear dispersion relation and dissipative spectral behavior. 
Models of similar complexity should be realizable with trapped ions. 
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I. INTRODUCTION 

Soft materials, such as granular media at vanishing external 
pressure or random polymer networks with loose connectivity, 
exhibit large deformations in response to applied mechanical 
perturbations. In both cases, there is a geometrical or topo- 
logical control parameter that determines the magnitude of 
the elastic moduli: the average overlap between grains or the 
mean coordination number of the network. By tuning these 
parameters one can reach a mechanical state characterized by 
a vanishing linear response, often termed sonic vacuum since 
the linear speed of sound vanishes^ When this happens even 
the tiniest strains propagate as supersonic solitons (or shocks) 
rather than ordinary phonons. The much studied (un)jamming 
transition is also marked by a progressive loss of rigidity that 
triggers a strongly non-linear response^ ^. At the transition 
these materials are not just soft they become fragile. 

The study of the impact of thermal fluctuations on these 
fragile states is still in its infancy. This is partly due to the 
fact that granular media and foams, that form the most natu- 
ral experimental arena to study these phenomena, are simply 
too large for thermal fluctuations to play any significant role. 
However, the existence of these athermal critical points rests 
on fine-tuning geometric parameters of their architecture that 
would be very sensitive to fluctuations if they were present. 

In this paper, we investigate a simple model that, at its 
critical point, displays a new quantum fragile behavior. As 
in the known soft counterparts, the mechanical response (the 
quadratic term in their Hamiltonian) vanishes, but in our 
model the mechanical properties are strongly afl'ected by zero 
point fluctuations in a way reminiscent of quantum phase 
transitions. The notion of fragility that we consider in the 
present stu dy is not necessarily related to the onset of a plas- 
tic response^^ES. Instead, we denote a classical or quantum 
mechanical system as fragile, if a parameter can be tuned 
such that the quadratic part of its Hamiltonian becomes sub- 
dominant so that its response to mechanical perturbations is 
intrinsically non-linear. Recent developments with trapped 
ions^ Rydberg atoms^^ and quantum optomechanics^^ 
give hope that synthesizing soft-matter inspired non-linear 
systems in the quantum regime may soon be within reach. 




FIG. 1 : (Color online) Sketch of the model Eq. ([TJ of a chain of quan- 
tum point particles in a local, pinning potential connected by springs. 
The blue, straight spring describes the harmonic term for the nearest- 
neighbor strain with strength A. The red, curved spring describes the 
Reggeon-type cubic nonlinearity for the strain with strength ig. See 
Ref. 9 for a feasibility study of a related Frenkel-Kontorova chain 
with trapped ions in optical potentials, Ref. 10 for small-polaron for- 
mation with trapped ions, and Ref. 11 for a trapped ion experiment 
of coupled quantized mechanical oscillators. 

The main physical result of our analysis is that the basic 
excitations at the critical point of a Reggeon ion chain (Fig.[T] 
Eq. ([T])), are neither ordinary phonons, that cannot exist be- 
cause of the vanishing rigidity, nor the strongly non-linear 
solitons that characterize similar critical points in the absence 
of fluctuations. Instead, quantum fluctuations give rise to me- 
chanical excitations with non-linear dispersion relations. This 
signature can be traced to the emergence of a length-scale de- 
pendent sound speed (or equivalently elastic modulus) that we 
study with the renormalization group to one-loop order. 

A. Model 

Specifically, we consider a non-unitary but PT -symmetric 
massive quantum chain with the Hamiltonian operator 

^1 A 

(1) 

describing N point-like quantum particles subject to a local 
potential of strength (arising for example in a deep com- 
mensurate pinning potential), that interact with their nearest- 
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neighbors with a harmonic-chain type coupling A and a cubic 
nonhnearity with (imaginary) coupHng constant ig (see Fig.[T] 
for a sketch). Here denotes the deviation from the equihb- 
rium position of the ^'th particle and is the conjugate mo- 
mentum operator fulfilling the discrete commutation relation: 
[p^,X£^] = -id a' (ft = I, lattice constant a=\). The Hamil- 
tonian in Eq. ([T]) is not Hermitian - it is not invariant under 
combined matrix transposition and complex conjugation. In 
addition to translational invariance in space i ^ i + n with 
n integer and time, Eq. ([T]) is invariant under a simultaneous 
parity transformation (P\ x -x) and complex conjugation 
or time-reversal (7~). 

Three noteworthy features of Eq. ([T]) are: 

(i) Reggeon-type, non-unitary field theories with cubic non- 
linearity have a successful history of describing geometric 
phenomena such as directed percolatioil^^^^ 

(ii) The cubic nonlinearity ig is chosen proportional to the 
strain similar to the granular chain ^ but with an integer- valued 
slightly above the Hertz law for spheres where the interaction 
scales as {x^ - x^T with a = 5/2?^ 

(iii) The initially non-zero mass from the pinning potential 
explicitly breaks the continuous shift invariance x^ ^ x^ -\- S 
thereby avoiding the Luttinger liquid fixed point of one di- 
mensional quantum systems^^.^^ 

The cubic nonlinearity preserves ^7~- symmetry and also 
simplifies the mathematical analysis, for example by giv- 
ing rise to anomalous dimensions already at the one-loop 
level. Bender and Boettcher^^ demonstrated that cer- 
tain, non-Hermitian, imaginary Hamiltonians, like the one 
in Eq. ([T]), can still describe sensible quantum mechanical 
ground states with real and positive spectra provided they 
fulfill a combined parity and time-reversal, ^7~-symmetry. 
Such !P7~- symmetric systems have since been realized ex- 
perimentally in optical waveguides^^ and it seems possible 
to engineer periodically structured nonlinearities as well^^. 
Related chains of oscillators with Hermitian nonlinearities 
such as the Fern ii-Pas ta-Ulam problenP^,^ quantum Frenkel- 
Kontorova chainJ^^Sl] chains with disordei^^, or chains under 
applied forces^ can also be used as model systems to access 
the strongly non-linear regime considered in the present study. 

We can recast the Hamiltonian operator in Eq. ([T]) as a La- 
grangian density. 
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FIG. 2: (Color online) Phases of the Reggeon quantum chain. The 
critical point marks the beginning of the fragile regime of small pin- 
ning potential (m^) and harmonic spring constant (A) relative to the 
nonlinearity {ig). 



B. Key results 

Our main motivation to study the somewhat artificial model 
Eq. ^ is that it displays some appealing critical features that 
may occur also in more complicated and realistic systems. At 
the same time, it is amendable to an approximate one-loop 
renormalization group (RG) treatment that yields the phase di- 
agram of Fig. [2] In particular, the RG can flow into the highly 
nonlinear regime close to the critical point, where violent fluc- 
tuations of the 0's lead to non-Luttinger liquid criticality for 
the mechanical excitations in the chain. It is well known that 
the presence of a continuous symmetry in combination with 
unitarity inevitably leads to Luttinger liquid behavior of con- 
ventional gapless quantum liquids in one dimension^^ Both 
of these constraints are relaxed in the highly non-linear regime 
of the non-unitary quantum chain. This is the reason why 
the imaginary cubic nonlinearity acts as a relevant perturba- 
tion to the Luttinger liquid behavior and leads to an unusual 
power-law scaling of correlators. At the critical point, when 
the chain becomes fragile, the mechanical excitations disperse 
as cl) ~ with z = 0.35 signaling a strongly momentum- 
dependent sound velocity. The associated spectral function 
becomes dissipative and displays a quantum critical contin- 
uum that may associated with the spontaneous breakdown of 
the ^7~- symmetry. 

The rest of the paper is structured as follows. In Sub- 
sec. 



II A we explain our RG set-up and show how the flow 



equations yield an interacting fixed point. The new physics 
associated with this fixed point is presented in detail in Sub- 
sec. InB] In Sec.|Ill| we draw our conclusions and point toward 



promising experimental systems. 



(2) 



that defines a lattice quantum field theory for a pseudoscalar 
living on a one-dimensional lattice with sites numerated by £ 
and T is the imaginary time coordinate ranging from to oo as 
appropriate for a zero-temperature quantum field theory. Fol- 
lowing Ref . |3T] we assume cp to transform as a pseudoscalar, 
that is, (p changes sign under space reflection then the in- 
teraction remains !P7~-invariant since / changes sign under T. 



II. ONE-LOOP RENORMALIZATION GROUP 

In order to perform our RG analysis we write the 
Lagrangian in a Fourier representation using 0^(t) = 
(pki^W^^^ ^ and analogously for frequencies oj. The lattice 
action in the thermodynamic limit is obtained as: 



Fao W= J ^ (p-k(-c^)[ + Aa, (1 - cosM) + ]cl)k{cD) 

8Ao 



^o r . 
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(3) 



where we have abbreviated f = f"^ ^ T ^ and f , , = 

dr J^^ ^ We have introduced already here the 

cutoff scale A e {Aq, 0}, where Aq is the "ultraviolet" mo- 
mentum scale of the order of the inverse lattice spacing at 
which the renormalization parameters Zaq , , and gA^ 
take their initial values. The endpoint of the RG flow is at 
A ^ and quantum fluctuations on all scales have been 
integrated into the efl'ective action T^lcp]- We exclude large 
momentum transfers from the analysis and approximate the 
trigonometric functions in Eq. ^ by their leading polynomi- 
als 1 - cos[^] ^ y -h ... and sin[^2 + ^3] ^ + ^3 + .... Now, 
the cubic vertex becomes proportional to the total momentum 
that flows through it. 

Note that the propagator in the quadratic term in the first 
line of Eq. ^ has a massive phonon dispersion, leading to 
spectral weight at 



(4) 



with Z, A, and attaining finite values. Throughout the mas- 
sive phase, the position of the peak will receive finite renor- 
malizations but otherwise the theory will still support well- 
defined quasi-particles. When approaching the critical point, 
for vanishing mass ^ 0, one immediately notices that a 
loop expansion of Eq. ^ will suff'er from logarithmic and 
power-law infrared singularities. We will see below that the 
RG flow to low frequencies and momenta in the presence of 
interactions gives rise to fractional gradients in space and time 
characterized by two difl'erent anomalous dimensions. 



A. RG Scheme 

To analyze the physics of these singularities, we employ 
the (formally exact) RG flow equation for the efl'ective action 
Fa [i//], the generating functional for one-particle irred ucible 
correlation functions in the form derived by WetterichP^ESl^ 
We mainly use this method due to its simplicity and versa- 
tility to tailor it to the problem at hand. In principle, similar 
results should be obtainable within an 6-expansion or a field- 
theoretic RG. We add a regulator Ra (specified below) to the 
quadratic part of our bare action Eq. ^ that introduces a cut- 
off' dependence into the eff'ective action so that T^lcp] smoothly 
interpolates between the bare action, Eq. ([3]), at the ultraviolet 
scale Fa^Aq [0] = Faq [0] and the fully renormalized eff'ective 
action in the limit of vanishing cutoff*: limA^o Ta [0] = F [0]. 
The Wetterich equation 



^ArA[0] = ^Trf- 



A^A 



2 'rfw^RA 



(5) 



has a one-loop structure and in a vertex expansion the /3- 
functions for the ^-point correlators are determined by (cut- 
off* derivatives of) one-particle irreducible one-loop diagrams 
with fully dressed propagators and vertices. The Tr is here just 
a frequency and momentum integration and F^^[0] = ^-^^ 
is the second functional derivative with respect to the fields 0. 



Za, Aa, 




FIG. 3: Diagrams for the flow equations. 



In the course of the RG ffow, new terms, not present in the bare 
action, will be generated. The objective of the present paper is 
to choose a simple truncation that reveals the qualitatively new 
physical features arising from the interplay of quantum ffuctu- 
ations and the imaginary nonlinearity. We therefore focus on 
the ffow of the scale-dependent phonon frequency renormal- 
ization factor Za, the momentum renormalization factor Aa, 
the mass and the non-linear coupling ^a- No local cubic 
interaction /gioc0^ is generated at one-loop. This is because 
the k-integration in the triangle diagram (see Fig.|3]) picks up 
a factor P and vanishes by symmetry. 

As an infrared regulator we use a mass-like cutoff* Ra = 
AaA^ so that ^A^A = 2AaA where, as customary, we drop 
the (higher-order) term ~ ^aAa in the scale-derivative of the 
cutoff*. Now, the infrared singularities at low frequencies and 
momenta of the phonon propagator 



GA(^,/:) = [Ff[0]+7?A 



Zaoj^ + Aa^^ -\-m\-\- Ra 



(6) 



are regulated without by-hand changing the "relativistic" scal- 
ing of the bare phonon dispersion oj ^ k. We now plug in 
Eq. d3|) supplemented by the mass and regulator terms into 
Eq. (|5]). By comparing coefficients of the ffelds in a vertex ex- 
pansion, we obtain the ffow equations for Za, Aa, m^, and gA- 
The corresponding diagrams are shown in Fig.|3] 

Upon using rescaled variables, the ffow equations can be 
brought into a simple form that facilitates comparison to the 
well known Wilson-Fisher ffxed point of the 0"^ theory below 
its upper critical dimension. The anomalous dimensions ab- 
sorb singular renormalizations to the frequency (Z) and mo- 
mentum (A) factor of the phonon propagator, respectively: 



rii = -^OaZa , ^ = -^OaAa . 

^A ^A 

Upon rescaling the mass and the cubic interaction as 



AaA2 



^A 



^A 



(J) 



(8) 
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we obtain their yS-functions: 
Ad^ml = {-2W^)fnl-lgl^^ (9) 

a^aIa = (-1 + 1< + + ^gi \ 2 • (10) 

^ ^ (^1 + thj^j 

The anomalous dimensions for the frequency and momentum 
factors, respectively, are different 



7 ^ ~9 



4 ^ ~2 



1 



1 



(11) 
(12) 



In the infrared A ^ 0, these flow equations have two fixed 
points. A trivial, non-interacting fixed point f a^o = with 
zero anomalous dimensions. More interesting is the infrared- 
stable, interacting fixed point ^ 0) with finite anomalous 
dimensions: 



jj^ = -0.04 , T/f = 0.61 . 



(13) 



We note that, as usual, the precise numerical values of the 
critical exponents depend on the renormalization scheme 
adopted, for example on the choice of the cutoff' function. 



B. New physics associated with interacting fixed point 



From the yS-functions Eq. ([9p2|), we see that the interact 



ing fixed point separates a massive phase from a region of the 
phase diagram (not studied in the present work and denoted by 
a dashed line in Fig. [2]) where presumably the "PT"- symmetry 
is spontaneously broken. 

The existence of a non-Gaussian fixed point with finite g* is 
ensured by two conditions: (i) the relevant dimensional run- 
ning leading to the -1 in the first bracket of Eq. (10) and (ii) 
the positivity of the second term ~ g^. In terms of power 
counting, (i) is actually similar to the conventional Wilson- 
Fisher fixed point in 3 dimensions, (ii) is a direct conse- 
quence of the imaginary nature of the cubic interaction in the 
Lagrangian Eq. ([2]). The same mechanism turns out to sta- 
bilize the real-valuedness of the fixed points of the local 
theory Another signature of the non-unitary interaction is 
a negative anomalous dimension as it was also pointed out in 
reference^ ^ and previously in the classic analysis of a 0^ the- 
ory by Fisher^^. We also find a negative t]^ for the frequency 
renormalization factor here in Eq. ^U) . 

A distinctive feature of the non-unitary quantum chain con- 
sidered here is that it is not relativistically invariant nor it 
seem to become so at the critical point, at least not within 
our one-loop approximation. The cubic vertex depends lin- 
early on spatial momenta (see below Eq. Q), since it couples 
neighboring sites. By contrast, it does not depend in the same 
way on frequencies because the action is completely local in 
time. While such frequency-dependent interaction will also 



be generated in the RG ffow, they appear only at higher-order 
in a loop expansion. As a consequence, at the one-loop or- 
der considered here, the anomalous dimension for the spatial 
momenta, rf^, is significantly larger and has the opposite sign 
to 7]^. At the critical point, the phonon propagator takes the 
form: 



1 



(14) 



resulting in a fractional value for the dynamical exponent at 
the interacting fixed point 



= 0.35. 



(15) 



This strongly modifies the phonon dispersion at the critical 
point to 



^critical 



rO.35 



(16) 



and leads to a critical continuum in contrast to the well- 
defined quasi-particle peak of the massive phase in Eq. (|4]). 



III. CONCLUSION 

In this paper, we presented a one-dimensional quantum sys- 
tem with a critical behavior that is diff'erent from a Luttinger 
liquid. To achieve this, we computed the one-loop renor- 
malization group ffow of a quantum harmonic chain supple- 
mented by a local mass term and an imaginary cubic non- 
linearity. The critical point we found may be viewed as a 
quantum fragile state of matter with unusual collective excita- 
tions and critical exponents. While we have no mathematical 
proof that the ^7~- symmetric icp^ quantum chain has a well 
deffned quantum mechanical ground state, it seems reasonable 
to extrapolate the work of Bender et al?^ on the local the- 
ory to a one dimensional lattice with translationally-invariant 
nearest-neighbor interactions. 

In the future, one could study local quantum quenche s of 
Eq. ([T]) in the strongly nonlinear regime to study solitons'^^^ 
or shocks^. We also work toward making synthetic quan- 
tum fragile matter accessible in future experiments by scal- 
ing up systems that have achieved nonlinear couplings at the 
single quantum lev el suc h as trapped ionP^, quan tum opto- 
mechanical system^l^'^' or optical waveguide 
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